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ABSTRACT 

We report measurements of geometric offsets between gas spiral arms and associated star forming 
regions in the grand-design spiral galaxy M51. These offsets are a suggested measure of the star 
formation timescale after the compression of gas at spiral arm entry. A surprising discrepancy, by 
an order of magnitude, has been reported in recent offset measurements in nearby spiral galaxies. 
Measurements using CO and Ha emission find large and ordered offsets in M5L On the contrary, 
small or non-ordered offsets have been found using the HI 21cm and 24/im emissions, possible evidence 
against gas flow through spiral arms, and thus against the conventional density-wave theory with a 
stationary spiral pattern. The goal of this paper is to understand the cause of this discrepancy. We 
investigate potential causes by repeating those previous measurements using equivalent data, methods, 
and parameters. We find offsets consistent with the previous measurements and conclude that the 
difference of gas tracers, i.e., HI versus CO, is the primary cause. The HI emission is contaminated 
significantly by the gas photo-dissociated by recently-formed stars and does not necessarily trace the 
compressed gas, the precursor of star formation. The HI gas and star forming regions coincide spatially 
and tend to show small offsets. We find mostly positive offsets with substantial scatter between CO 
and Ha, suggesting gas flow through spiral arms (i.e., density-wave) though the spiral pattern may 
not necessarily be stationary. 



1. INTRODUCTION 

The density wave theory (|Lin fc Shulll964| ) has been a 
central para digm for the fo rmation of spiral structures 
in galaxies. iRobertsI (|1969| ) suggested that star forma- 
tion (SF) in spiral galaxies is triggered by a spiral den- 
sity wave and predicted that Ha emission around newly- 
formed stars should be offset from and appear down- 
stream of the gas spiral arm due to flow through the 
density wave. The offsets would not only be evidence of 
the density wave, but would allow measurement of the 
angular speed o f the spiral pattern and the timescale for 
star formation ([Egusa et al.ll2004[) . given a velocity dif- 
ference between the matter and the spiral pattern and 
assuming that the SF timescale is constant. The offset 
and circular velocity can be measured observationally as 
a function of radius, and should obey the relationship, 

A9(r)^ B - (fi(r) - n p )At (1) 

where A6(r) is the azimuthal offset angle between gas 
spiral arms (A) and star forming regions (B), fi(r) is the 
angular velocity of the gas, ft p is the angular pattern 
speed, and At is the time that the gas takes to evolve 
into young massive stars after spiral arm entry. 

A discrepancy has emerged in recent stud i es of offsets 
(ITamburro et all 120081 lEgusa et all 120091: IFovle et al l 
1201 lh . In an analysis of ~ 10 galaxies. lEgusa et alJ(|2009L 
hereafter E09) derive d offsets larger by an order of magni- 
tude in general than ITamburro et al.1 1)20081 T08). More 
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surprisingly. IFovle et al.l (|2011l Fll) reported primarily 
non-ordered offsets compared to those predicted by the 
standard density- wave theory (Eq. 1). This discrepancy 
directly affects the estimates for il p and At, and thereby 
their physical implications for star formation. Indeed, 
the measured star formation timescales are inconsistent: 
5-30 Myr for E09 versus 1-4 Myr for T08. The shorter 
timescale may imply very short lifetimes for giant molec- 
ular clouds (GMCs), i.e., their exceedingly rapid forma- 
tion and destruction and associated rapid star formation 
(T08), while the longer timescale may be consistent with 
the gravitational collapse of gas at typical giant molecu- 
lar cloud gas densities (E09). 

This study explores the potential causes of the discrep- 
ancy. Two major differences among these studies are 
the choices of emission tracers and t he methods o f offse t 
measurement. E09, and orig inally lEgusa et al.1 (|2004[ ). 
used CO (J=l-0) and Ha as tracers of the dense gas 
and recently-formed young stars, respectively. T08 and 
Fll used HI 21 cm emission and IR 24/im emission to 
trace the locations of the compressed gas and young mas- 
sive stars, respectively. T08 and Fll also compared low- 
resolution CO and 24/im emission on a subset of galaxies. 
Fll also investigated other tracers such as UV emission 
to trace SF and 3.6 /an emission to trace the old stellar 
population. CO is the established tracer of the dense star 
forming molecular gas, while HI, if concentrated around 
spiral arms, may coincide with the dense gas. Ha ob- 
servations typically provide a higher angular resolution 
than observations at 24/im. However, Ha observations 
could miss dust-obscured star forming regions, whereas 
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24/im emission suffers very little from extinction. 

The methods of offset m easurement could a lso be the 
cause of the discrepancy. lEgusa et al.l (|2004D and E09 
compared the location of emission peaks by eye. This 
"Peak Tracing" Method is intuitive and easy to apply, 
but is possibly biased. T08 and Fll used an automated 
cross-correlation method to avoid any human bias. This 
"Cross- Correlation" Method is apparently less subjec- 
tive, but surprisingly its application leads to inconsistent 
results even when the same data are used (T08 and Fll). 

The offset measurements may constrain emerging the- 
ories of spiral arm formation. In iV-body simulations, 
the spiral arms appea r to form as mat erial arms in- 
stead of density waves (jWada et al.l 120 111) . Spiral struc- 
tures are constantl y formed and broken apart by grav- 
itational shearing (Goldrcic h fc Lvnden-BelH Il965| ). In 
this case, no clear offset would be expected across the 
spiral arm. The stellar and gas spiral arms are most 
likely co-spatial, the gas being pulled toward local stel- 
lar potential minima, instead of passi ng through them a s 
predicted by the density - wave theory ()Wada et al.ll2011| ). 
iDobbs fc Pringlel (|2010T l performed hydrodynamic simu- 
lations in stationary and transient stellar spiral struc- 
tures (i.e., density- wave spiral versus material spiral, re- 
spectively) and found that the offset depends strongly on 
the nature of the st ellar spiral st r uctur es. The ordered 
offsets measured by lEgusa et~al] (j2004ft . T08, and E09 
support the density wave theory in its simplest form with 
a single pattern speed, while the result of Fll could indi- 
cate a transient spiral structure without a fixed pattern 
speed. Resolving the discrepancies in the offset measure- 
ments is therefore increasingly important. 

In this paper, we reexamine the measurements of 
the geometric offsets by T08, E09, and Fll, using the 
previously-adopted four tracers, HI, CO, 24/mi and Ha 
emission, and using both the Peak Tracing and Cross- 
Correlation methods. We focus on M51 to clarify the 
primary causes of the discrepant results from previous 
studies. We conclude that the discrepancy comes from 
the different tracers, mainly because the HI 21 cm emis- 
sion traces the gas dissociated by young stars rather than 
the parental gas for star formation. We will demonstrate 
that the two measurement methods provide consistent re- 
sults if appropriately applied. In addition, we will discuss 
some caveats in applying the offset model for deriving the 
pattern speed and the SF timescale. 

M51 is among the galaxies analyzed by T08, E09 and 
Fll. However, they did not find consistent geometrical 
offsets AO (and, almost equivalently, the star formation 
timescale At). For example, T08 measured At = 3.4 ± 
0.8 Myr (fi p = 21 ± 4kms" 1 kpc" 1 ), while E09 obtained 
At = 13.8±0.7Myr (fi p = 40 ± 4kms _1 kpc" 1 ), i.e. 
a factor of 4 difference in At (a factor of 2 in A 
summary of the tracers used and results can be found in 
Table [U 

The data for the gas and star formation tracers are 
described in §2. The offset measurement methods of the 
previous studies are discussed in §3. Our results with 
four emission tracers, HI, CO, 24/im, and Ha, are dis- 
cussed in §4. The measurement methods and tracers are 
compared in §4.2 and §4.3, respectively. The cause of 
the discrepancies among previous measurements and the 
limitations of the model are discussed in §5. A summary 



of this work is given in §6. 

2. DATA 

In order to elucidate the causes of the discrepancy, we 
compare the tracers (i.e., HI, CO, 24/im and Ha) pre- 
viously used in T08, E09 and Fll The HI data are 
from the HI Nearby Galaxy Survey (jWalter et al J 120081 
THINGS), an HI survey made with the Very Large Array 
at the National Radio Astronomy Observatory. We use 
the data reduced with natural we ighting. The CO data 
arc from iKoda et al.l (pool 12011 and were obtained by 
combining Combined Array for Research in Millimeter 
Astronomy (CARMA) and Nobeyama 45m telescope ob- 
servations as a part of the CARMA-Nobeyama nearby 
galaxies (CANON) CO (J=l-0) survey. The 24/xm and 
Ha dat a are from the Spitzer Infrared Nearby Galaxies 
Survey (|Kennicutt et al.ll2003l SINGS). The 24/zm obser- 
vations were made with the Spitzer Space Telescope and 
the Ha observations were taken with the 2.1m telescope 
at The Kitt Peak National Observatory. The spatial res- 
olutions of HI, CO, 24/im and Ha data are ~6" (280 pc 
at 9.6 Mpc), 4" (160 pc), 5.7" (256 pc), and 1.9" (88 pc), 
respectively. 

The HI and 24/im data are the same as used in T08 
and Fll. The CO data have a higher resolution and sen- 
sitivity than the data fro m the BIMA Survey of Nearby 
Galaxies (BIMA-SONG; iHelfer et al.l l2003t i that were 
used in T08 and E09; however, both resolve the molec- 
ular spiral arms in M51, and there should not be much 
difference in the offset analysis. Fll used much lower res- 
olution (13") CO (J=2-l) data fro m the HERA CO-L ine 
Extragalactic Survey fHERACLES; iLerov et ail2009ft in 
their analysis. The Ha image is not the same as the 
one that E09 used, though it has a similar resolution. 
We do not expect a significant difference between the 
two, because the analysis is weighed significantly to high 
signal-to-noise regions in bright spiral arms. A slight 
improvement in the S/N does not change the positions 
of the bright peaks. We adopt the following parameters 
for M51: a distance (D) of 9.6 Mpc, a major axis posi- 
tion angle (P.A) of 22° and an inclination angle (i) of 
20° (jSofue et al.ill999f ). We will test the different sets of 
these parameters and confirm that the choice does not 
affect our conclusions. 

3. OFFSET MEASUREMENT AND ANALYSIS 

The procedure for the offset method has two major 
steps: (1) measurements of azimuthal offsets between gas 
and star formation tracers, and (2) the determination of 
a star formation timescale and pattern speed using the 
measured offsets and a rotation curve. Our primary fo- 
cus is to investigate possible causes of the discrepancy 
in step (1). The discrepancy between the previous offset 
measurements (i.e., large offsets in E09, small offsets in 
T08, and no ordered offsets in Fll) could be due to dif- 
ferences in the methods employed and/or in the adopted 
tracers of dense gas and star forming regions. We briefly 
summarize the two methods here. 

Both methods are bas ed on theoreti cal predictions of 
the density- wave theory (|Robertslll969l ). Assuming a flat 
galactic rotation curve and a constant pattern speed, the 
relative velocity between the gas/stars and the spiral pat- 
tern changes with radius. If the timescale of SF after spi- 
ral arm compression is a constant, this differential mo- 
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tion results in a spatial offset between the gas spiral arm 
and star forming regions, and the offset changes with ra- 
dius. If we measure the offset as a function of radius, 
the SF timescale a nd pattern speed can be determined 
(jEgusa et al.l 120041 ) . Figure Q] shows two model plots of 
the expected offset versus galactic radius, assuming SF 
timescales of 1 and 10 Myr and other parameters appro- 
priate for M51 (i.e., pattern speed and rotation velocity). 
Figure [T] (a) shows the azimuthal offset in the galactic 
disk, while Figure Q] (b) shows the offset in the sky. The 
azimuthal offset is expected to be large at small radii and 
approximately zero at large radii. 

For both methods, the four images, HI, CO, 24/im 
and Ha, are regrided to the same pixel size us- 
ing the Multichannel Image Reconstruction, Image 
Analysis, and Display (MIR IAD) software package 
(jSault. Teuben fc Wrightlll995| ) and are deprojected us- 
ing the MIRIAD task "deproject" . The deprojected im- 
ages were transformed into phase diagrams (radius, az- 
imuth) using the MIRIAD task "polargrid". After the 
production of the phase diagrams, the analysis continues 
differently for each method as discussed below. 

We wrote our own IDL programs to reproduce the dif- 
ferent measurement techniques of E09, T08, and Fll, 
following the detailed procedures discussed in these pa- 
pers. 

3.1. Peak Tracing Method 

A measure ment method f or the geometrical offset was 
developed by lEgusa et"al~l (|2004D and updated in E09. 
The offset is measured by eye as the azimuthal angular 
offset between the emission peaks of gas and star forming 
region tracers. This Peak Tracing method starts by bin- 
ning the phase diagrams in the radial direction with an 
appropriate bin size to match the spatial resolution of all 
tracer images (here we use a bin size of 5" ) . The emission 
peaks are visually determined in an azimuthal intensity 
profile of each bin. The offsets are defined as the az- 
imuthal shifts of the peaks between the two tracers. We 
find that this visual identification tends to pick up only 
peaks with a high signal-to-noise ratio. Obviously, some 
radii are not included in the analysis if no emission peaks 
could be identified in their azimuthal profiles. Figure [5] 
shows examples of identified peaks in all the four tracers 
in a part of Arm 1: (a) HI , (b) CO, (c) 24/rni, and (d) 
Ha. Gas is flowing in a counter-clockwise direction. 

3.2. Cross- Correlation Method 

T08 cross-correlates the azimuthal profiles of the gas 
and SF tracers in each radial bin. They define a normal- 
ized cross-correlated function, 



^k[{xk - x)(y k -i - y)} 
y/'Ekixk - x) 2 Y, k (y k - y) 2 



(2) 



where x% and y k are the fluxes in gas and SF trac- 
ers (in T08 and Fll, HI and 24 /im, respectively), x 
and y are the mean values of the fluxes. I is the lag 
(or the shift) in pixel units between the two profiles. 
The denominator normalizes the cross-correlation func- 
tion based on the variance of each data set; this normal- 
ized cross-correlation function is useful when adopting a 
cutoff value to remove low signal-to-noise data at differ- 
ent wavelengths. In principle, the cross-correlation func- 



tion shows a maximum at the value of I = ^ max that best 
aligns the two profiles - finding the peak of the cross- 
correlation function should be enough to obtain Z max . In 
practice, we fit a polynomial to the profile of the cross- 
correlation function around the peak to reduce the effect 
of noise. T08 uses a fourth-degree polynomial, and we 
follow their method. Examples of the cross-correlation 
function for a radial bin (at the galactic radius of 2.51 
kpc) are shown in Figure [3] for each set of tracer com- 



binations. The offset is simply AO(r) = Z r 



T08 



adopted a lower threshold (i.e., 0.2) for the normalized 
cross-correlation function to reject measurements with 
low significance and Fll used a threshold of 0.3. We set 
the threshold for an acceptable cross-correlation value to 
0.3 to avoid low signal-to-noise measurements. 

3.3. Parameters and Errors 

We apply both the Peak Tracing Method and the 
Cross-Correlation Method to two gas tracers (CO and 
HI) and two SF tracers (Ha and 24/im). The phase dia- 
grams are binned radially in 5" increments for all data. 
This bin size is chosen to match the resolutions of the 
HI and 24/im data. We denote the spatial offset be- 
tween the gas (tracer A) and star forming regions (B) as 
AOa^b and measure AQ(r) HI ^ 2 4fim, AQ(r) HI ^ Ha , 
AO(r)co^24/itm, and A8(r) co ^ij a . If multiple peaks 
exist in the cross correlation function, the largest is cho- 
sen, since by definition this is the location where the 
peaks of the spiral arms are best aligned. For the Peak 
Tracing Method, the locations of any multiple peaks are 
checked against maps of M51 to make sure that they cor- 
respond to t he spiral arms inste ad of emission in the inter 
arm regions (jKoda et al.l 120091 ) . 

Two types of error could be introduced in the anal- 
yses: (1) due to uncertainties in determination of peak 
positions in each tracer, and (2) due to the misassoci- 
ation of the peaks in gas and SF (e.g., in the case of 
multiple peaks). The first type of error is easy to esti- 
mate and is simply related to the resolution of the im- 
ages in both methods. This relation may not be so ob- 
vious in the case of the Cross- Correlation Method, but, 
for example, if there are point sources convolved with 
point-spread functions (PSF; i.e., resolutions) the cross- 
correlation function is simply a multiplication of the two 
PSFs with a variable lag. Its width (or the associated 
error) is therefore determined by the sizes of the reso- 
lutions. In all the data, we measure only the positions 
of bright peaks and the positional accuracy should be a 
small fraction of the resolutions. We here conservatively 
adopt half of the resolution as the error. This choice pro- 
vides an error similar to that estimated by E09, but much 
larger than that of T08. The second type of error is dif- 
ficult to estimate, and we neglect it for now. This error, 
however, should appear as scatter in our measurements. 

4. RESULTS 

We qualitatively compare the images of the different 
tracers before making quantitative measurements. To 
distinguish the two arms, we define Arm 1 as the spiral 
arm that is further from the companion galaxy and Arm 
2 as the one that connects to the companion. 

Figure H] shows phase diagrams of the four tracers. 
Contours of star formation tracers (Ha and 24/zm) are 
shown on gas tracer images (CO and HI). Gas flows in the 
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direction of increasing azimuthal angle. Figure [5] shows 
the relative locations of the peaks in the four tracers with 
respect to CO (i.e., green line). If we compare the peak 
locations of the SF tracers, Ha and 24/im, with respect 
to CO peaks, the Ha emission appears downstream of 
the CO emission when an offset is seen. This is consis- 
tent with the findings in E09. The 24/zm peaks show a 
similar trend and appear to be leading the CO peaks. 
Most important are the locations of the HI peaks as they 
are mostly at the leading side of CO, the well-established 
tracer of star-forming gas. The HI peaks appear to be 
spatially closer to the peaks of the two star formation 
tracers. This trend is inconsistent with the hypothesis in 
T08 and Fll that the HI traces the compressed, parental 
gas for SF. The Ha and 24/xm peaks appear at similar 
locations in most radial bins, although occasionally they 
show offsets. The offsets between CO and Ha/24/zm are 
evident over a range of radii, while HI shows little or no 
offsets from Ha and 24/im emission. Another notable 
point is that CO is predominantly in the inner part of 
the galaxy, while HI is in the outer part. 

We note that the spiral arms do not form narrow 
smooth lines; on closer inspection they are fragmented 
and wiggle back and forth between the leadi ng and trail- 
ing si des (most obviously seen in CO; see also lKoda et al.l 
2009). This causes an additional intrinsic error in the off- 
set measurements. 

We apply the Peak Tracing Method and Cross- 
Correlation Method to these four data sets, and quantita- 
tively verify the cause of the discrepancies in the previous 
measurements. We omit the very central part (radii < 2 
kpc), since the spiral structure is tightly wound and we 
cannot reliably associate the star-forming gas and resul- 
tant young stars in that region. 

4.1. Comparisons to Previous Works 

Figure [3] compares the offsets measured using the same 
tracers and methods as the previous studies. The left 
column shows the results from the Peak Tracing method 
applied to one spiral arm (i.e., Arm 1; as in E09). The 
right column shows the same offsets but from the Cross- 
Correlation Method applied to both arms (as done in T08 
and Fll). The offsets between HI and 24/im are shown 
in the top panels, and the offsets between CO and Ha 
are shown in the bottom panels. 

The measured offsets in Figure [5] are consistent with 
the previous measurements of T08, E09, and Fll in 
terms of the ranges of their amplitudes. E09 compared 
CO and Ha data using the Peak Tracing Method and 
found offsets up to 25-30 ° in Arm 1. Figure [5] (bottom 
left) shows our corresponding measurement. The ampli- 
tude range of our offsets between CO and Ha is similar 
to that in E09 (see their Figure 5, Arm 1). Our range 
in Arm 2 is also consistant with that of E09 (see Fig- 
ure [5]) . E09 found irregular offsets between CO and Ha 
along Arm 2, and our measurements for Arm 2 see simi- 
lar irregularities, which are discussed later in this section. 
T08 analyzed HI and 24/im using the Cross- Correlation 
Method (see their Figure 4) and found smaller offsets 
(mostly less than 5-10 °), except for the innermost radii 
(some large offsets of 15-20° < 2 kpc). In our compari- 
son of HI and 24/im (Figure [SJ top right) we find similar 
small offsets, most of which are under 10 ° in > 2 kpc, 
consistent with T08 results. We were unable to measure 



offsets at the innermost radii (discussed above). Fll also 
measured offsets between HI and 24/im using the Cross- 
Correlation Method and found that they are mostly un- 
der 10° (their Figure 10), which is comparable to our 
amplitude range. The amplitude of the offsets found by 
T08, E09, and Fll agrees with our measured offsets when 
we use the same set of tracers and measurements method. 
Our measurements reproduce the amplitude discrepancy 
in the previous studies. 

Figure E] shows plots of all of the azimuthal offsets mea- 
sured for each set of tracers as a function of galactic ra- 
dius for the Peak Tracing method (left two columns) and 
the Cross-Correlation method (middle two columns). We 
measure the offsets in the two arms separately in these 
four columns. T08 and Fll analyzed the two arms to- 
gether to measure one offset at each radius. For compari- 
son, we also analyze the two arms together and show the 
plots in Figure |6] (right column). E09 found a strange 
trend in Arm 2, where the offset increases with radius 
(as opposed to the prediction in Figure [1]) . Our results 
for Arm 2 from the Peak Tracing method show similar 
behavior. Along Arm 2, the emission from the SF tracers 
appear to be scattered over a large area (though mostly 
at the leading side of the arm). Therefore, we discuss 
only Arm 1 in the following analysis when we analyze 
the two arms separately. 

T08 also analyzed CO and 24/im data (their Figure 7) 
and found that the offsets are even smaller than those 
seen between HI and 24/im. We found the contrary. The 
offsets between CO and 24/im in Figure [6] ( third row) are 
larger compared to those between the HI and 24/im (top 
row). Fll also compared a lower-resolution CO image 
with a 24/xm image and noted possible evidence for or- 
dered offsets in M51, though they did not quantify the 
CO— >24/im offsets further. Indeed, compared to their 
Figure 11, we find even larger offsets at the small radii 
(Figure [S](£op right)). Even if we smooth our CO data to 
a 13" resolution to match Fll's, we still find larger offsets 
than Fll. Fll also analyzed the two arms separately, but 
report no difference. Again, this is inconsistent with our 
results. Part of the reason for this difference may be that 
they used higher transition CO (J=2-l) data while we 
used CO (J=T0). The higher transition em ission might 
be ex cited due to heating by young stars (jKoda et al.l 
[2012]) . If that is the case, CO (2-1) should appear closer 
to SF tracers (i.e. smaller offsets). 

4.2. Method Comparison 

The two measurement methods provide, in general, 
consistent results and do not appear to be the cause of 
the discrepancy. Figure [7] directly compares the mea- 
sured offsets in Arm 1 from the Peak Tracing Method and 
the Cross-Correlation Method (i.e., the points in Figure 
[6]) and shows good agreement, with some outliers. The 
dotted lines have a slope of one. Two solid lines in each 
panel show ±5° from the dotted line. We find that 70% 
of the data points are within ±5°. Except for some out- 
liers, most of the measured offsets are comparable, and 
neither method appears to bias the offsets in any sys- 
tematic way in any of the panels. The two methods give 
roughly consistent offsets. The top left panel only has 
3 out of 19 radii where the two measurement methods 
do not show a correlation. For the top right panel 6 out 
of 19 points do not agree. The bottom left panel has 5 
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points out of 19 and the bottom right panel has 6 out of 
17 points that do not show a correlation. 

E09 discussed that the two arms show quantitatively 
different offsets and analyzed them separately with the 
Peak Tracing Method. The difference between the two 
arms is clear in Figure [6] (first and second columns). 
Since the two methods generally provide consistent re- 
sults when applied to single arms, it is natural to expect 
differences in the measured offsets when the two arms 
are analyzed simultaneously using the Cross-Correlation 
Method. We see this in our analyses of both arms simul- 
taneously (Figure ®(fifth column) and separately (Figure 
[B] (third and forth columns)) with the Cross-Correlation 
Method. We note that both T08 and Fll analyzed 
the arms simultaneously, while Fll also isolated indi- 
vidual arms; the difference between the two arms was 
not mentioned in either case. From our work, regard- 
less of the measurement method (the Peak Tracing or 
Cross-Correlation Methods) and the analysis (the two 
arms separately or simultaneously), the offsets are mostly 
positive, but do not appear ordered. 

4.3. Tracer Comparison 

The choice of tracers appears to be the dominant cause 
of the discrepancy. Large differences are seen among the 
panels in Figure [5] ( left column) and H] (middle column), 
where the two arms are analyzed separately. Most re- 
markably, HI and 24/im - the combination employed by 
T08 and Fll show the smallest, if non-zero offsets for 
M51. Similarly, very small and zero offsets are seen be- 
tween HI — > Ha. If we use CO as a tracer of the dense, 
compr essed gas (e.g., CO — > Ha used by lEgusa et al.l 
(2004) and E09), the offsets are much larger. Clearly, 
the difference in gas tracer is the dominant cause. 

Figure [8] shows a direct comparison of the offsets 
between gas tracers and SF tracers in a given radial 
bin. The left column shows the offsets from the Cross- 
Correlation Method, while the right column shows those 
from the Peak Tracing method. The top panels compare 
the offsets found using HI with Ha and CO with Ha. 
The bottom panels show the same, but for 24/im. The 
offsets between HI and the SF tracers are mainly zero; 
however at the same radii the offsets between CO and the 
SF tracers are larger. In addition, when non-zero offsets 
are measured for both gas tracers, the offsets between 
CO and the SF tracers appear larger than the offsets be- 
tween HI and the SF tracers. This is evident in Figure 
[8] as most data points appear in the top-left parts of the 
panels. Occasional outliers exist, but there are no more 
than two in each panel. 

There are multiple components that could be con- 
tributing to the HI emission, which would give differ- 
ent offsets (and possibly non zero offsets) depending on 
which component is dominant at the radius in question. 
For example, HI may trace some of the gas compressed 
upon entry into a spiral arm, resulting in positive offsets, 
but would also trace the gas photo-dissociated by newly- 
formed young stars, leading to zero offsets. There is also 
an extended background component in the HI emission 
that is virtually everywhere across the galactic disk. We 
can also find some spatial displacements between HI and 
CO peaks in Figure 2]- HI is mostly downstream of CO, 
indicating that the HI emission is not simply from the 
dense gas, but multiple sources. This trend is also seen 



in Figure (SJ where AOco^hi is plotted against radius. 
The positive offsets suggest that HI is downstream of the 
CO. We will discuss this more in £15-11 

As for SF tracers, 24/im and Ha peaks are mostly con- 
sistent; therefore, the offsets with both tracers are more 
or less the same in Figured! However, there are some no- 
ticeable differences at some radii. For example, in Figure 
[4] there are radii at which Ha peaks appear downstream 
of 24/im peaks. This is perhaps due to the fact that the 
distribution of 24/im emission is affected by the distri- 
butions of both dense gas/dust and illuminating stars. 
If newly-formed stars escape from the dense parental gas 
(as we see in Figure 2]), the 24/im emission may not trace 
the locations of the young stars any more, because the 
peak of dust density is left behind. The relevance of the 
tracers to the offset analysis is discussed in 35.11 

4.4. Miscellaneous Differences 
4.4.1. Resolution 

The offset measurement could be sensitive to the reso- 
lution of the data. From Figure [TJ we can determine that 
the azimuthal offset is expected to be 5 — 30" in angu- 
lar distance on the sky (1.5 — 6kpc), if a star formation 
timescalc of 10 Myr (E09) is assumed. This range is mea- 
surable with the spatial resolutions of the data. On the 
other hand, if the timescalc is 1 Myr as suggested by T08, 
the spatial resolution required to see offsets would be 
1 — 5" . Therefore, the resolutions of some of the adopted 
data would be marginal (i.e., 4" in CO, 6" in HI, 1" in 
Ha, and 6" in 24/im) if At is as small as found by T08. 

In order to test for a possible resolution dependence, 
we smooth the Ha data to match the resolution of the 
24^m data (~6") and re-measure the offsets, and find 
practically no appreciable difference from the measure- 
ments at a higher resolution for HI— >Ha and CO— >Ha. 
A visual comparison of the Ha smoo thed maps with the 
24/im maps shows that the two emissions are very well 
aligned and trace similar emission sources. The differ- 
ence in resolution of the SF tracers does not appear to 
be effecting the measured offset. 

We also smooth the CO data to the resolution of the HI 
data (6") and re- measure the offsets with the SF tracers. 
Again, we find no significant difference in these measure- 
ments with the higher resolution images for CO— )-24/im 
and CO— >-Ha. Therefore, the resolution dependence is 
negligible for the data that we analyzed. 

4.4.2. Bin Size 

We also test two different radial bin sizes (5" and 2" 
bins) in measuring offsets. These are the bin sizes used 
in T08 and E09, respectively. We find little difference 
between the offset patterns for all four combinations of 
the tracers; and hence, this is not the cause of the dis- 
crepancy. 

5. DISCUSSION 

5.1. Nature of Gas and Star Formation Tracers 

We found that the most significant cause of the dis- 
crepancy between the previous studies is the difference 
in the gas and star formation tracers used. We also found 
that CO emission is the best tracer of the compressed gas 
in spiral arms and that Ha is the best tracer of the as- 
sociated star forming regions. Thus, they are the most 
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useful for the offset measurement among currently avail- 
able archival data. The previous studies adopted differ- 
ent tracers for tenable and practical reasons. Here we 
discuss the advantages and disadvantages of each tracer. 

Ha emission is from HII regions around young massive 
stars whose lifetimes are short (< lOMyr). There is no 
doubt that discrete Ha peaks pinpoint the locations of 
very recent star forming regions. In addition, many sets 
of Ha data are readily available in archives, and typically 
have a higher resolution than the Spitzer 24/im data, an- 
other widely-used tracer of star formation. On the other 
hand, Ha emission could be easily obscured, especially in 
dense star-forming gas in spiral arms. We might be pref- 
erentially finding HII regions far away from the parental 
gas, biasing the measured offsets toward the large side 
(T08). This, however, has proved not to be the case. 
The ability of Ha emission to trace the locations of star 
forming regions has been shown recently b y a compar- 
ison of Ha and Paschen a images of M51 ()Egusa et al.l 
1 2 1 If) . Even in the densest regions in a spiral arm, the 
majority of HII regions, if not all, can be seen in Ha 
emission, though their fluxes could be significantly at- 
tenuated. Ha emission is the best locater of HII regions 
for offset measurement since their positions are the only 
parameter required. 

The 24/xm emission is less sensitive to extinction and 
might be a better locator of star forming regions that 
are embedded deeply in a dense star-forming spiral arm. 
However, it may not be as straightforward a star for- 
mation tracer as previously thought, since the emission 
could be fro m dust heated by longer-lived, older stellar 
populations (|Liu et al.ll20lH ) . In addition, the dust could 
be heated by the continuum emission longward of the 
Lyman limit from young stars. The shorter-wavelength 
emission is mostly absorbed in the vicinity of the young 
stars and forms HII regions, but the longer- wavelength 
photons can potentially travel farther and heat up dust 
not immediately associated with the young stars. If this 
is the case, the 24/im may not pinpoint sites of recent 
star formation. The 24/xm image is also sensitive to the 
distributions of both heating sources and the gas/dust. 
If the gas spiral arm and recently-formed stellar spiral 
arm are spatially offset (e.g., the stellar arm should b e 
downstream within a corotation radius; lRobertsl fl969), 
the young stars may illuminate the gas/dust spiral arm 
from the front-side and shift the peak of the 24/im emis- 
sion. There is little doubt that young stars are the domi- 
nant source of dust heating around star forming regions, 
and their discrete appearance is very clear in the 24/im 
image. However, contamination by the escaped photons 
from HII regions and the older stellar population could 
be possibly biasing the offset measurements using 24/xm 
emission to the small side. 

A concentration of CO (J = 1 — 0) emission is the 
clearest tracer of dense molecular spiral arms and is 
aligne d perfectly with n arrow dust lanes in optical im- 
ages (|Koda et al.l 12009). The critical density for col- 
lisional excitation of the low- J CO line is high (e.g., 
- a few x 100 cm" 3 iScoville fc Sanders! Il987l) . and all 
young star forming regions in the Milky Way are asso- 
ciated with giant molecular clouds (GMCs), which are 
bright in CO. There is no doubt that high-resolution CO 
data show the locations of the gas spiral arms that are 



being compressed by the spiral density- wave. 

It also had seemed reasonable to assume that the HI 
emission is adequate to locate the enhancements of gas 
density in a spiral arm (T08 and Fll). However, the 
small, and often zero, offsets of SF tracers from HI, as 
opposed to the large offsets from CO, are a clear sign 
that the HI emission peaks are not tracing the com- 
pressed gas in the spiral arm. In fact, the HI peaks are 
almost always at the downstream side of the CO peaks 
(Figure [2] and Figure HJ) . The zero offsets to SF tracers 
most likely indicate that the HI emission is traci ng gas 
that is photo-disso ciated by recent star formation (jAllenl 
120021) . Blitz et al.l (|2007j ) showed that GMCs are often 
found at the peaks of HI emission in the Hi-dominated 
galaxy M33, but there is also more extended HI emis- 
sion across the galaxy in regions without GMCs. HI 
peaks may trace locations where the gas will not form 
into stars. It is remarkable that the fraction of molec- 
ular gas does not change a zimuthally across spiral arms 
in M51 (jKoda et al.l 120091) . and therefore the effect of 
photodissociation is small in terms of the global evolu- 
tion of the gas phase in the spiral galaxy. However, this 
small change makes a huge impact in the identification 
of the HI arm - since overall the abundance of atomic 
gas is low - and likely caused the discrepancy in the 
previous studies. Perhaps there is some enhancement of 
HI emission along the gas spiral arm, but in most cases 
it is washed out by the stronger emission of the photo- 
dissociated HI g indicated by the zero/small offsets 
between HI and SF (T08 and Fll). These observational 
results require a reconsideration of the notion of a short 
GMC life-cycle, an argument largely based on the off- 
set measurements with HI and 24/im emissions (T08). 
The offset measurement does not indicate the lifetime of 
GMCs, since the gas stays mostly molecular even after 
spiral arm passages and star formation. 

Therefore, we conclude that CO and Ha are the best 
tracers of the compressed gas in spiral arms and the 
associated star forming regions, respectively. The HI 
emission does not necessarily trace the compressed, star- 
forming gas in the spiral arm; instead it traces predomi- 
nantly the gas photo-dissociated by recent SF. Emission 
from 24/im could be tracing an older stellar population 
in addition to areas of recent SF and many not pinpoint 
the site of SF due to the offsets between the dust and 
young star distributions. 

In this work our discussion of the offset is limited 
to M51, and we can not make any general conclu- 
sion. Nevertheless, our finding that HI traces the photo- 
dissociated gas, rather than star-forming dense gas, offers 
a natural explanation for the general discrepancies in the 
previous studies. The best combination of tracers (CO 
and Ha) shows clear spatial offsets and indicates larger, 
mostly positive offsets than previously suggested using 
the HI data (T08), although the positive offsets are not 
ordered as predicted by the standard density-wave the- 
ory. 

5.2. Model Limitations to the Offset Measurement 

The offset method is based on only a few assumptions, 
i e., a constant f2 p and At and pure circular rotation. Its 
simplicity assures robustness, but has some weaknesses 
as well. We note and summarize several limitations here 
which are likely related to the large scatter found in the 
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plots presented above. Some of the limitations come from 
the model (Equation 1), and the others are from intrinsic 
location-to-location variations of gas conditions and star 
formation in spiral arms. 

An intrinsic spatial variation of CO and Ha distribu- 
tions along spiral arms is a source of a systematic error. 
As discussed in gl the spiral arms are not simple, contin- 
uous structures; at small scales they appear as ensembles 
of more discrete clumps in both CO and Ha images. The 
spiral arms go back and forth between the downstream 
and upstream sides, and this wiggling directly affects the 
offset measurement. The typical amplitude of the wiggle 
is about 200 pc and contributes to the scatter in Figure 
[6l We should note that this type of error tends to can- 
cel out and does not bias offset values systematically. In 
addition, the molecular gas shows filamentary structures 
(or spurs) in the i nterarm regions a t the downstream 
side of spiral arms (|Koda et all 12009). These interarm 
structures can potentially limit the offset measurement, 
especially using the cross-correlation method where such 
individual structures are neither identified nor rejected 
in the analysis. 

One of the other limitations is the finite lifetime of 
HII regions. The star formation timescale At derived by 
the previous studies ranges between < 3Myr (T08) to 
> lOMyr (E09). 

Some of these timescales are as short as (or less than) 
the typical lifetime of HII regions (~ lOMyr). The de- 
rived timescales are meaningful only when they are com- 
parable with or longer than the the lifetime of HII re- 
gions. If one obtains a shorter timescale, more careful 
consideration, such as taking into account the age dis- 
tribution of HII regions, is necessary before any physical 
interpretation of the results. 

The assumption of a constant pattern speed is com- 
monly adopted, but it is possible that this assumption is 
not valid. Indeed , it has been sugge sted that f2 p changes 
radially in M51 (|Meidt et al.ll2008D . In particular, the 
grand-design spiral arms in M51 c ould be dr i ven b y the 
companion galaxy (NGC 5195). lOh et all (|2008l ) dis- 
cussed, based on theoretical modeling, that Q p in tidally- 
driven spiral arms may vary with radius. The mostly 
positive offsets that we found in M51 are an indication 
that the material is flowing through the spiral pattern. 
Therefore, the spiral needs to be a wave, not a concen- 
tration of material This, however, does not mean that 
the spiral pattern speed is constant with radius (i.e., the 
traditional density wave) . If the stellar spiral density en- 
hancement lives longer than the gas crossing time, the 
observed results can be explained without a stationary 
spiral pattern. Note that the assumption of constant f2 p 
has an impact neither on the offset measurement itself 
nor on our conclusion that the cause of the discrepancy 
in the previous studies is the choice of tracers. 

The circular rotation of galaxies is another assump- 
tion that most studies adopt. T08 calculated the direc- 
tional changes in gas orbits due to the spiral arm poten- 
tial and concluded that its impact on At is very small 
(i.e., substantially less than a factor of two). Another 
important factor is the change of the amplitude, not 
only the direction, of gas velocity due to a symmetric 
galactic potential. Ga s orbits in spiral galaxies are of- 
ten described as ova ls (|Wadalll99"llOnodera et al.1120041; 
iKoda fc Sofudl2006l) . In an axisymmetric galactic poten- 



tial the gas moves inward and outward; the speed being a 
maximum at the points on the semi-minor axis of the or- 
bit, and a minimum on the semi-major axis. Spiral arms 
form primarily due to this slowdown on the semi-major 
axis; the gas stays longer on the spiral arms, which leads 
to the density enhancement. More precisely, the spiral 
density enhancement occurs even if we do not account 
for the shock (which produces even more enhancement), 
and is proportional to the passage times through the spi- 
ral structure t cro ss °c l/(f2(r) — tt p ). Q(r), and there- 
fore t cross , change substantially during the orbital mo- 
tion, and can change the s urface density in spira l arms 
by an order of magnitude (|Onodera et all 12004 their 
figure 11). 

To estimate the effect of this slow down to At, we 
estimate a difference in rotational velocities between cir- 
cular orbits assumed in Eq. (1) and elliptical orbits that 
form spiral arms. Assuming that the gas is traveling in 
an elliptical orbit with the eccentricity e = yl — b 2 /a 2 
(maximum and minimum distances a and b, respectively) 
and that the galactic potential is isothermal (i.e., flat 
rotation curve, $ oc ln(r)), conservation of energy and 
angular momentum provide 

v a /TT 2 , (e 2 - 1) . > 

where v a is the rotation velocity at the maximum dis- 
tance a for an elliptical orbit, v c is the circular rota- 
tion velocity at radius a. f2 a and O c are the angular 
speeds corresponding to v a and v c . For example, as- 
suming an elongated orbit (e = 0.6-0.8), derived pat- 
tern speed (fip = 30kms _1 kpc -1 ), and rotation veloc- 
ity (~ 200 km s -1 ), At oc l/(ft — £l p ) changes from its 
value under circular rotation by a factor of 2-3 (but up 
to 10-30) at radii of 4-6 kpc . If the non-circular mo- 
tion, the very essence of spiral arms, is taken into ac- 
count, the derived At under the assumption of circular 
rotation is an underestimate by a factor of 2-3. It is 
important to note that the non-circular motions do not 
explain the discrepancy in the measured offsets. They 
affect the interpretation of the azimuthal offsets, but not 
their measurements. [Note again that the main focus of 
this study is the measurements.] In order to confirm our 
simple analytic calculations, we compare the simplistic 
cir cular rotation mo d el wit h the numerical simulations 
bv IDobbs fc Pringld (|2010D and Fll that include non- 
circular motion. We use the offsets measured by Fll be- 
tween the gas and 100 Myr old clusters in the simulations. 
The pattern speed of their stationary spiral model is 20 
kms -1 kpc -1 and the rotation curve is roughly flat with 
a peak velocity of 220 kms -1 . We adopt these parame- 
ters for our circular rotation model. Figure [TU] shows the 
expected offsets for the circular rotation model, Eq. (1), 
compared to the offset from the numerical simulations. 
It is clear that the circular rotation model overpredicts 
the offsets by a factor of two in this radius range. There- 
fore, the gas motion across the spiral arm slows down in 
the more realistic numerical simulations. The difference 
is expected to be larger around the co-rotation radius (~ 
11 kpc in this parameter set). These results are consis- 
tent with our simple calculation of non-circular orbits. 
We conclude that the circular rotation model adopted 
in the previous studies tends to underestimate the SF 
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timescale by a factor of 2-3. Development of an orbit 
model that fits M51 is beyond the scope of this paper, 
but this discussion provides a caveat that the At derived 
from the offset method is likely a lower limit of the SF 
timescale, which could be about an order of magnitude 
greater when the non-circular motion is taken into ac- 
count. This could also be evidence against a short GMC 
lifetime. 



6. CONCLUSION 

The discrepancy found between the previous offset 
measurements (T08, E09, Ff f) is significant. E09 found 
relatively large offsets, which are consistent with the pre- 
dictions of the density-wave theory and possibly gravita- 
tional collapse of the dense gas formed in spiral arms. 
T08 found offsets an order of magnitude smaller, which 
would suggest extremely rapid star formation and de- 
struction of giant molecular clouds in a spiral arms (i.e., 
faster than the free-fall timescale of the gas). Most strik- 
ing are the non-ordered offsets found by Fll, which could 
argue against the traditional density-wave theory. To 
elucidate the cause of the discrepancy, we repeated the 
previous studies using the gas tracer emissions (CO and 
HI) and SF tracer emissions (Ha and 24/im) and apply- 
ing the two measurement methods that were developed 
in the previous studies. 

We analyzed the spiral arms in M51 and found that 
the primary cause of the discrepancy is the use of dif- 
ferent gas tracers. In particular, the HI 21 cm line 
emission traces predominantly the gas photo-dissociated 
by recently-formed stars, but not necessarily the com- 
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We compared several other differences in the analyses 
of the previous studies, and found that nothing but the 
tracers contribute significantly to the discrepancy. The 
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measurement methods, spatial resolution, and bin size. 
The dominant cause of the discrepancy is the different 
choice of the tracers in the previous studies. HI emis- 
sion does not necessarily trace the sites of compressed 
gas nor star formation, and contaminated significantly 
by photodissocation due to young stars. 
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TABLE 1 

Summary of Previous Offset Studies 
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Fig. 1. — The expected offsets determined from Eq. (1) for a constant v ga s of 200 kms -1 and C p of 30kms —1 kpc -1 . Two star formation 
timescales are assumed, Ai = 1 Myr and 10 Myr. (a) The azimuthal offset in the galactic disk plane as a function of radius. (6) The angular 
offset in the sky. The offset is expected to be largest at smaller radii and approaches zero with increasing radius. 




Fig. 2. — Images of a portion of Arm 1 in (a) HI, (b) CO, (c) 24^m, and (d) Ha. The positions of peaks identified by the PT method in 
radial bins (centered on the concentric arcs) are plotted on the images. The black circles in the lower-left corners indicate the resolutions 
of the data. The gas flow direction is counter-clockwise. As a reference, the CO peaks are connected with green lines in all panels. 
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Fig. 3. — Examples of the cross-correlation function for the four combinations of the tracers at a radius of 2.51 kpc. The dashed curve 
show a forth degree polynomial fit around the maximum peak. The dashed blue line shows the location of the maximum peak. The offset 
value is listed in the top right of each panel for each tracer combination. 
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Fig. 4. — Phase diagrams of combinations of the four tracers for comparison of our results with those in T08, Fll (top-right) and E09 
(bottom- left). The offsets are measured using the methods adopted in T08, Fll (top panels) and E09 (bottom panels). 
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Fig. 5. — Azimuthal offsets between gas and SF tracers. The top row shows the offsets between HI and 24fim emission (A^hi— >24/im)> 
while the bottom row shows those between CO and Ha (Ailco^Ha)- The Peak- Tracing (PT) method is applied to Arm 1 for the left 
column, and the Cross-Correlation (CC) method is applied to both spiral arms simultaneously for the right column. The top-right panel is 
directly comparable to the measurements by T08 and Ff 1 (as it uses the same tracers and method). The bottom-left panel is comparable 
to the measurement by E09. 
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Fig. 6. — Azimuthal offsets between gas and SF tracers as a function of radius. The two spiral arms are analyzed both separately and 
simultaneously with the Peak Tracing (PT) method and Cross-Correlation (CC) method. Left two columns: the offsets determined by the 
PT method with the two arms analyzed separately, Middle two columns: by the CC method with the two analyzed separately, and Right 
column: by the CC Method with the two analyzed simultaneously. The rows from top to bottom show the offsets between HI — > 24/im, HI 
— ¥ Ho, CO — > 24/im, and CO — > Ho. The model predictions of offsets for star formation timescales of 1 and 10 Myr are overplotted (see 
Figure [TJ. 
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Fig. 7. — Comparisons between the offsets using the Cross-Correlation (CC) and the Peak Tracing (PT) methods. Offsets are measured 
at each radial bin for Arm 1, (see Figure |B). For all four tracers, the offsets measured with the two methods are in general agreement, 
with some outliers. Neither method appears to be biased in any systematic way. The fraction of data that show correlation is ~ 70%. For 
reference, a dashed line with a slope of one is plotted and two solid lines are plotted to show ±5° . 



CC Method 



PT Method 





30 








20 


8 


X 




t 


' 


O 




O 





© 






- 1 n 



-4 


/ 

^ / 

\ 


/ 

" / 







30 
20 
' 




AO HI ^Ha (°) 



4 


B 

jji^tt ; 


/ 

■ / 







-10 10 20 30 
AG HI ^24^,m (°) 



A0 HI ->Ha (°) 



-10 10 20 30 
AO HI ^24,am (°) 



Fig. 8. — Comparisons between the offsets measured using different gas tracers on Arm 1. The top row shows the comparisons between 
Af2co->Ha and Af2jji— >Hm while the bottom row shows those between ACcO— >24/^m and Af2ni— >24/^m- The CC method is applied in the 
left column, and the PT method is applied in the right column. For reference, the dashed line has a slope of one (i.e., the two offsets are 
equal) and the two solid lines show zero (no offsets). Most points appear in the top-left of each panel, indicating that the offsets between 
CO and SF tracers are larger than those between HI and SF. The offsets between HI and SF tracers are often zero. 
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Fig. 9. — Offsets between CO and HI as a function of radius. Most offsets arc positive, indicating that the HI peaks are on the downstream 
side of the CO. 
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Fig. 10. — Comparison of offsets from two models: the analytical circular ro tation model adop t ed in most offset studies (dash-dotted line) 
and more realistic numerical simulations using a stationary spiral pattern by Dobbs & Pringlc (2010) and Fll (crosses). The simulations 
include non-circular motions across spiral arms, and therefore trace any slowdown of gas and stellar motions. The two models adopt the 
same global pattern speed, Q p = 20 kms" 1 kpc -1 , and a flat rotation curve with a peak velocity of 220 kms -1 . The offsets are measured 
between the gas and a 100 Myr-old stellar population. The offsets in the numerical simulations are lower by a factor of ~ 2 than those 
predicted by the circular rotation model, indicating that the gas and stars slow down considerably during spiral arm passage. 



